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1. INTRODUCTION 

Recently, various techniques and method have been developed and proposed for the numerical 
solution of ordinary differential equations and partial differential equations. The hyperbolic Goursat PDE arises 
in linear and nonlinear equations involving mixed derivatives. Researchers embark on developing various 
numerical methods to solve the problem such as Runge Kutta method [1], Two Dimensional Differential 
Transform Method [2], Finite Difference Method [3, 4], Variational Iteration Method [5] Adomian’s 
Decomposition Method [6] and Explicit Numerical Scheme [7]. 

Fuzzy transform has been utilized and introduced as an approximation method in order to solve 
various differential equation including ordinary differential equation and partial differential equation [8-15]. It 
has been demonstrated that if the functions are replaced by an approximation scheme, at some point where 
complex computation appear, it can be simplified by fuzzy approximation method. 

The main idea of this paper is to apply the fuzzy transform at both sides of Goursat hyperbolic PDE 
and approximate the functions using finite difference method and arithmetic mean. 


2. PRELIMINARIES 
This section is divided into two subsections, which are fuzzy transform and hyperbolic Goursat partial 
differential equation. 
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2.1. The fuzzy transform 

Fuzzy transform has been introduced by Irina Perfiliev [16]. The Fuzzy Set Theory concept by Lotfi 
Zadeh [17] is the fundamental of development of fuzzy transform. Some basic definitions is included to 
demonstrate the Fuzzy Transform. 


b-a 


Definition 1 Let x, =a+h(i-1) be nodes on [a,b] where h = i 
n — 





n= 2 and i =l,...,n. We 


say that fuzzy sets A,(X),...,A,(x) defined on [a,b] are basic functions and form a fuzzy partition of [a,b] 


if the following conditions hold true for eachi = 1,...,n: 


A, :[a,b] > [0,1], A, (x,) =1, 


A, (x) =0 if XE (eysta) where Xq = a Xn = b, 
A(x) , : 
is continuous, 
A, (x) Legs, | [X; Xi] 


is strictly increases on and strictly decrease on 


Di 4) = a alt X45] 

A(x, -x) =A,(x, FD sera xe [0,A] i=2,...n—-Ln> 2 

A,,,(x) = A(x—h) pra xelath,b], i=2,...,.n-2,n>2, 

Definition 2 Let AD) yie.s A, (x) be basic functions which form a fuzzy partition on [a,b]. The n-tuple 


of real numbers [F, sete F | given by 


[POA Cdr 
f= i 
[ A,(x)dx 


=1,....n 


is the direct fuzzy transform of function f with respect to A,(x),...,A, (x) . The [Feck®, | are called the 
component of Fuzzy Transform. 


Definition 3 Let the [F, yoees Fe | be the direct Fuzzy Transform of function f with respect to 


A,(x),...,A, (x) . Hence, the function 


f= DRAG). 


is called the inverse Fuzzy Transform of f. 


Definition 4 Let f(x, y) be an arbitrary continuous function on D=[a,b]x[c,d] and real matrix 
Ppt 2 
F°TfI=([Fy lm given by 


LI LEMACOB Oddy 
[Pawson 


y 





is the Fuzzy Transform of f with respect to the given fuzzy partition. The reals F, are the components of the 


Fuzzy Transform of f. In the next section, a brief explanation on the hyperbolic Goursat partial differential 
equation is discussed. 
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2.2. Hyperbolic goursat partial differential equation 
The Goursat equation is a second order hyperbolic PDE involving two independent variables x and 
y . The general form of the Goursat problem is given by Wazwaz [18]. 


6*u(x, y) 
Oxy 

u(x,0) = (x), 

u(0, y)=y(y), 

$(0) = y (0) =u(0,0) 

O<x<a, O<y<b 


= g(x, y,U,uU,,U,), 


2 
Ou(x,y) , deers 
where —————— iis the mixed derivative in space x and y while the term a(x, y,U,U, ,U, )is a function 
X 


consists of a combination of the independent variables x and y , the dependent variable u and the derivative 


terms u, and uy. 


The Goursat problem arises in linear hyperbolic PDE and nonlinear hyperbolic PDE. A linear Goursat 
PDE necessarily satisfies two particular conditions otherwise, it will be considered as a nonlinear equation. 


The first condition is the power of the dependent variable uw and the derivative term wu, and WU, must be one 


[18]. Besides, the coefficients of the dependent variable 1 and the derivative term uw and U y must be the 


constant or the independent variables x and y [18]. In this paper, a Goursat PDE will be experimented. Some 


source of errors might occur during handling the Goursat PDE. The round off error usually occurs while 
allocating with the computing software because of its inability to deal with certain numbers. The terms must 
be approximated by using any certain suitable method where the truncation error might take place as well. 
Every approximation at certain grid points will be computed by using a computing software. The propagation 
error might occur since the approximation from present iteration depends on the approximation from the 
previous iteration and so forth. Hence, several difficulties in discretizing the Goursat equation with the 
derivative term may arise. Thus, it is aneed to propose a new method on hyperbolic Goursat partial differential 
equation. This new method is discussed in the next section. 


3. THE PROPOSED OF FUZZY TRANSFORM METHOD FOR HYPERBOLIC GOURSAT 
PARTIAL DIFFERENTIAL EQUATION 

In this section, a new numerical method involving fuzzy transform is developed. Certain 
simplification of complex computations could be achieved by replacing the original function with an 
approximation model [8]. The main idea is to apply the fuzzy transform to both sides of the hyperbolic Goursat 
PDE. The application of fuzzy transform to the hyperbolic Goursat PDE lead to the transformation from 
differential equation to the algebraic equation which is solvable. The aid from MATLAB is very beneficial in 
calculating the exact solution and approximate solution at each grid point as well as the relative errors and the 
average. The relative error is used to describe the accuracy of the numerical method. The proposed numerical 
method can be extended in solving differential equation involving more variables. 


Let the domain D? be given by a Cartesian product of two real intervals X =[O,l]and Y =[0, R] 


Let u(x, y) be a continuous solution of the following hyperbolic Goursat partial differential equation. 


O°u(x, y) _ 


B(x, y,U,u,,U,) 
OxOy * 


After applying the fuzzy transform, the hyperbolic Goursat PDE is transferred to the following 
algebraic equation: 
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F*[u"]=F’[g] 














Where: 
Un Uy Ui, 
Ua, Mag, So Loe 
F?[u my = 7 7 ‘ 
Oe, ee ee 
O°u(x, y) 
is the matrix of the fuzzy transform components of —————— , and components of g represented by matrix 
Xt 
of the fuzzy transform as below: 
Gj, G). - G,, 
Gy Gy Gy, 
2 . . . 
Eel= 
IG, cn oc G,,, | 
eR: 5 OUOE I). ee te} 
ince the partial differential equation —————— is unknown, the partial derivative is replace e 
S the partial differential equat A kn the partial d ti placed by th 
Xt 


forward finite difference approximation. 


Oru(x, y) . 
OxOy 


is replaced by 
—j lua th’, yth?)—u(xt+h", y)—u(x,y +h") +u(e, y)] 


Next, the approximation of U;; is as follows: 


PP a.008 oraray 
U» = 





[/ [CoB Qydxay 





Mite pr lMr thy thule hy) uy HR) + UC yi}aons, (y)dxdy 


~ 
~ 


[ [CoB (dray 
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1 [ [ACB Odedy [ [A @B, Qddxay 


is wa. yth’ A, (x)B ;(y)dxdy . i [uc y))A,(x)B, (y)dxdy 
[/ [AB Q)dxdy [/ [CoB ()dxdy 

















1 





G+D G+) ny, ~ Uicisy = U,,) 


> 


In order to produce more accurate approximation of G,., an arithmetic mean technique has been 


L 


utilized. The approximation is as below: 


8(X, y,U,U,,U,) is replaced by 
| ie : x 7 
qle@+h yth’)+e(xy)t+g(xth y+ g(x,y th] 


d eb 
[ J, eG. y.u.u,.u, A, (OB, (ydedy 


. d pb 
| [ A, (x) B , (y)dxdy 





ij 


1 ; ; 
PE (Glec iy sis cone acer y+ eauy +h 9], 008, (nara 


~ 
~~ 





[{[ 4. @oB, Cardy 
[Pec rny +h CoB onardy [fe ACB, Orardy 
[FAB aray [{ [A CoB, Oday 
Eilean my CoB ondrdy PII (eos y+), CB, ddedy 
[J ACB, O)axay [' [4 @oB, Oddy 











1 
= 7 Gosnien + Gi + Goan + Cs) 


Thus, the recursive equation is as below 


a 
nh? 


1 
( wan ~F aayy —Giguy + Yi, ) “4 (are + G+ Gos + Gigs) 


4. NUMERICAL EXPERIMENT 
Consider the following hyperbolic Goursat PDE 


U,, =U, 
u(x,0) =e", 
u(0, y) =e’, 


O<x<4, O<y<4. 
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The exact analytical solution of the Goursat problem is u(x : y) =e*'” [19]. The accuracy of the 


method is investigated by observing the relative error and average relative error at different values of step size, 
has shown in Table 1, Table 2 and Table 3. 


Table 1. Approximate Solutions at Various Grid Points and Step Size, h 








, : Step size, h 
U(xy) Bee solution 0.0025 0.005 0.025 0.05 
U(0.25,0.25) 16487213 1.6487214 16487216 1,6487298 16487554 
U(0.50,0.50) 2.7182818 2.7182823 2.7182837 2.7183280 2.7184667 
U(0.75,0.75) 4.4816891 4.4816905 4.4816949 4.4818339 4.4822684 
U(1.00,1.00) 7.3890561 7,3890598 73890708 73894227 73905227 
U(.25,1.25) 12.1824940 12.1825023 12.1825272 12.1833242 12.1858159 
U(1.50,1.50) 20.0855369 20.0855545 20.0856072 20.0872951 20.0925718 
U(.75,1.75) 33.1154520 33.1154876 33.1155944 33.1190128 33.1296999 
U(2.00,2.00) 54.5981500 545982199 54.5984295 54.6051377 54.6261111 





Table 2. Relative errors, E i,j (h) with h=0.0025 


x 








y 0.5 1.0 1.5 2.0 
0.5 1.6996495e-07 2.8515340e-07 3.6294990e-07 4.1532779e-07 
1.0 2.8515340e-07 4.9606933e-07 6.5061317e-07 7.6293317e-07 
15 3.6294990e-07 6.5061317¢-07 8.7527807e-07 1.0485725¢-06 
2.0 4.1532779¢-07 7.6293317e-07 1.0485725¢-06 1.2796936¢-06 





Table 3. Average relative errors, E(h) with several step size, h 





Step size, h Average relative error, E (h) 
0.0025 4.2688467e-07 
0.005 1.7111414e-06 
0.025 4.3506679e-05 
0.05 1.7771664¢e-04 





It indicate that the relative errors become smaller and approaching zero as the step size, h tends to 
zero. Hence, the accuracy of the method point out to be improved as the step size, h decreasing. For a better 


understanding, the average relative errors, E(h) for the step size h=0.0025, h=0.005, h=0.025 and h=0.05 are 
shown in Table 3. Figure | shows the graph of exact solution and approximate solution. 


Exact Solution Approximate Solution at h=0.005 


Z-axis 











Figure 1. The graph of exact solution and approximate solution at h=0.005 
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The graph of approximate solution at h=0.005 is graphically similar to the graph of exact solution. 
Thus, we deduce that the relative errors of approximate solution are significantly small. Hence, both 
presentations for the approximate solution either in a tabular form or graphical illustration proved that the 
method is an effective method. 

The fuzzy transform method is successfully applied to hyperbolic Goursat PDE where the numerical 
solution is accurate by comparing the relative error at different step size. The significant of fuzzy transform is 
it reduce the complexity of differential equation by transforming it to the algebraic equation. The proposed 
numerical method can be extend to differential equation involving more variables. The advantage of the 
proposed method is for every term in algebraic equation, it can be approximated by different existing 
approximation techniques and successively developed a new numerical method. Thus, the fuzzy transform 
method is a good numerical method for solving hyperbolic Goursat PDE. 


5. CONCLUSION 

In this paper, a new method has been successfully developed by utilizing the Fuzzy Transform for 
solving various engineering problem represent by hyperbolic Goursat partial differential equation. The fuzzy 
transform is utilized on both side of Goursat equation to produce a continuous approximation. Since various 
engineering problems involve complex model, the complexity can be reduced by applying fuzzy mathematics 
concept. The difficulties and complications of the proposed numerical method is overcome by developing a 
MATLAB program as well as with the aid of MAPLE software. Therefore, the computational numerical 
solution has been presented in this paper. Numerical methods are very significant to solve complex model 
considering the numerical methods can be effectively demonstrate in cases dealing with complex problems for 
which analytical solutions cannot be obtained or hand calculations cannot be made. Hence, this method can be 
an alternative numerical method to solve any physical problem represented by a hyperbolic Goursat PDE. 
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